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Abstract 

We numerically study the interacting quantum Hall skyrmion system based 
on the Chern-Simons action. By noticing that the action is invariant un- 
der global spin rotations in the spin space with respect to the magnetic field 
direction, we obtain the low-energy effective action for a many skyrmion sys- 
tem. Performing extensive molecular dynamics simulations, we establish the 
thermodynamic phase diagram for a many skyrmion system. 

PACS numbers: 75.10.-b, 73.20.Dx, 64.60.Cn 
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The application of a strong magnetic field to a two dimensional electron gas "quenches" 
the kinetic energy of the electrons - removing it from consideration by constraining them to 
cyclotron motion. This allows other physical effects such as interactions and correlations to 
dominate. |1J When the number of electrons, N e , is an integer multiple of the number of flux 
quanta, N^, the groundstate of the system is a set of completely filled Landau levels. If we 
define the filling factor v = N e /N^, then for v = 1, 1/3 and 1/5 the groundstate is fully spin 
polarized. This polarization arises not simply from the Zeeman splitting but rather from 
interactions; indeed it persists even in the limit of zero Zeeman coupling. f^-Q The system 
minimizes its Coulomb energy by minimizing the overlap of electrons in the many-body 
wavefunction. When the electrons all have the same spin, the amplitude for two electrons 
to be at the same spot must vanish due to the exclusion principle, thereby reducing this 
electrostatic energy. 

Near v = 1 the low energy excitations of the system are distortions of the spin orientation, 
which are coupled to the local charge density. The system arranges spins of adjacent electrons 
so that they are nearly parallel, forming smooth distortions or textures in the spin field. 
This can be described by a non-linear sigma model, which predicts localized charged 
excitations consisting of dimples in the spin field called skyrmions. It follows 

that if one moves away from v = 1 the groundstate will include such textures as (possibly 
pinned) quasiparticles. This has been verified experimentally by NMR measurements of the 
spin polarization. p~0j , pTTJ| Other experiments such as optical absorption |12|] and transport 
13fl are consistent with this picture. 



At higher densities (i.e. further away from v = 1) these quasiparticles themselves start 
to interact with each other. Each skyrmion has an overall U(l) symmetry corresponding 



to rotation about the z-axis. [T^J] When two are brought near each other, the lowest energy 
configuration is to align themselves anti-ferromagnetically. When many skyrmions interact, 
there is a competition between the Coulomb interaction, which favors a triangular lattice, 
and this antiferromagnetic XY interaction, which favors a bipartite one such as a square 
lattice. The resulting phase diagram is quite rich, and has been investigated theoretically. 



While the skyrmion lattice has not been seen directly, some experiments are suggestive. 
NMR measurements display a dramatic increase in the nuclear spin relaxation rate as a 
function of | v — 1| , |jIT[ and the low temperature heat capacity of multilayer 2DEG samples 
has a dramatic peak that may be related to the melting of a skyrmion lattice. [Tj| This 
anomalous, large (xlO 5 ) peak in the specific heat in GaAs quantum wells is postulated to 
come from a coupling of all of the nuclear magnetic moments of all the atoms adjacent to 
the quantum well, due to the symmetry breaking field of the Skyrme crystal. 

In this paper we numerically study the phase diagram of the Skyrme crystal. Elsewhere 



we have derived a low-energy effective action, for the many body skyrmion system. [19 



This effective action accurately describes the renormalization due to the higher Landau 
level mixing and is shown to be accurate for large and small skyrmion sizes. Here we use 
this action to variationally determine a two-body interaction potential between skyrmions. 
We then use this two-body interaction in a molecular dynamics (MD) simulation method 
to study the thermodynamic phase diagram of a many body interacting skyrmion system. 
Unlike earlier studies, we do not presume an a priori lattice structure. We incorporate the 
full, long range nature of the interaction allowing for competition between incompatible 
lattice structures and the possibility of frustration. At zero temperature, we find that the 
triangular lattice is formed below V\ = 1.025 and the square lattice becomes stable above v\. 
At high filling factor, the square Skyrme crystal will melt due to strong quantum fluctuations. 



| 14| Finite temperature simulation shows that at v = 1.1, the square Skyrme crystal melts at 
about T m = 0.035p s . For the appropriate choice of experimental parameter, i.e. p s = 1.5K, 
T m is estimated to be 50mK, which is very close to the value given by experiment. [|TI| 
We begin with the effective action Se for many skyrmion system |1 

^[m] = i£^)|J 1 2 + £^|J S r 



+ jjdr (Vm) 2 + ^Bra 2 
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+ ^A {0) -J*--E_L^- k xJ s (1) 

v , v f k 

where k = p/m*, V e s(k) = V(k)/(1 + bk) is a screened Coulomb interaction between 
skyrmions, V(k) = 2ne 2 /ek is a bare interaction, e is the dielectric constant, b = £ 2 /cib 
with a B = (h 2 e/m*e 2 ), i = (hc/\e\B)^ 2 is the magnetic length, and p s is the spin stiffness 
e 2 /16^/2ne£. The skyrmion three-vector J* is related to m(r), the local orientation of the 
spin texture via: 

J t = g^W(<^ m x ^m) • m . (2) 

By noticing that the presence of a finite Zeeman gap makes the skyrmion action invariant 
under global spin rotations by an angle <ft m the XY-spin space [0, we chose to describe a 
skyrmion as a point particle with charge q s = ve and phase <fi. The Lagrangian for a many 
skyrmion system can then be simplified 

+ ^A ( °»-J s -7E^-kxJ s . (3) 

The skyrmion transport mass M s and the moment of inertia I s are given by the following 
relations 

M s = — [ drf\JZ\ 2 (4) 
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2nm*£ 2 r u 

l^dr — (m x a y m±) ■ (5) 



where S represents an area occupied by a single skyrmion and a y is a Pauli spin matrix. 
The function V(R, x) is the interaction between two skyrmions, where R = |Rj — Rj| is the 
relative distance and x = (0i — <fij) relative spin orientations. The last two terms in Eq.(|3]) 
take into account the exchange statistics of a skyrmion. [Tj| For small Xi the two-body 



interaction can be decomposed into the following suggestive form 



V[R, x\ = V (R) + V A {R)[l + cos(x)]. (6) 

At large distances, Vq(R) goes like cfe/eR as a bare Coulomb repulsion between two point 
charges and Va(R) decays exponentially due to a finite Zeeman gap. Note that the angle 
dependent interaction Va{R) favors an antiferromagnetic spin alignment, i.e. \ = tt, trying 
to form a bipartite lattice, while Vq(R) prefers a triangular lattice which does not belong 
to bipartite lattices. This two-body interaction is difficult to calculate variationally: A 
configuration of two skyrmions must be allowed to relax without allowing the skyrmions to 
migrate apart. We chose an ansatz that fixed the relative angle between two skyrmions 
to be x an d the distance to be R, and then calculated the minimum energy configurations 
numerically. Our variational form is: 

2 e i0(x)/2 



u{z,z*) = K2/ 



z 2 



X 2 



e -\z-R/2\/2£ _|_ e -\z+R/2\/2£ (^) 

where z = x + iy denotes the position vector, £ = (2-k/E z ) 1 I 2 with E z = (l/2)gfiB / p s , 
and A is a variational parameter which controls the extent to which the spin fields are 
deformed. The spin fields m(r) are related to the complex field u(z,z*) via the following 
relation: m x + im y = 2a>/(l + |c<;| 2 ). |7[] When the distance R is much larger than £, our 
variational ansatz describes two nearly free n=l skyrmions, and near z = R/2, u(z,z*) = 
\(z - R/2)/\ s ]e lz ~ R/2l/2 t with X s = X 2 /(2R) and 9(x) chosen to be 0. The variational 
energy of an n=l skyrmion based on the above ansatz yields 1.789 in the units of 47rp s 
and A s = 2.34£, which is remarkably close to the exact value 1.788 from the field theoretic 



action. [19| In the opposite limit of R = 0, it depicts an n=2 skyrmion, whose variational 
energy is calculated to be 3.939, which is very close to the exact value of 3.913 in the units of 
4np s , which lends a strong credence to our variational ansatz. By minimizing the action in 
Eq.(|I]) for the proposed two skyrmion ansatz with 9{x) = 0, the two-body interaction Vq(R) 
is calculated. In order to calculate the angle-dependent interaction Va(R), we adjust A to 
minimize the energy for the following ansatz for the gradual spin twist 9(x), which rotates 
m with respect to the z axis: 

9{x) = 9 t&nh(2x/R) 



which uniformly inserts twist between the two skyrmions. Here ir — 9q is the relative phase 
difference \ between two skyrmions. This ansatz does not change the Zeeman energy, since 
m z stays the same. In Fig. (TTJ) we plot the interaction strengths V(R,tt — 9$) as a function 
of distance R and relative angle 9$ between the textures. We can break this into two parts 
as described in Eq.(^) and the interaction profile thus obtained is given by 

Tr/ „, „ 0.35 + 0.175.R 

^" 3 ' 578+ 1 + 0.0213^ + 0.0147^ (8) 

V A (R) = 0.Q29e-° AUR (9) 

where Vo and Va are in the units of 47rp s , R in the unit of t. One can notice that Vq(R) 
bends over to a certain finite value as R decreases due to the soft core nature of skyrmion 
and d 2 Vo/dR 2 becomes negative at v c = 1.2. This leads to a complex pattern at filling 



factors around and above v c [20|, which will most likely be preempted by quantum melting. 



TM Based on the interactions obtained above, we employ overdamped molecular dynamics 



simulations at finite temperature 



3& 



- J2 ^Vk(Ry) [1 - OOs(0y - It)] + £?(t) 

V+^T = ~ E Vk(B^) sin(^ - vr) + g(t) (10) 

mmm)= 2 v«T5 ij 5(t-t'). (n) 

Here rj a is the damping parameter with a — R, <f), = Rj — R^, 0^ = 4>i — 4>j, and £f (t) is 
a Gaussian white noise which represents the thermal fluctuations following the fluctuation- 
dissipation theorem. It can also help the system to reach the true ground state by providing 
a random noise. 

The number of skyrmions iV was fixed to be 576 and the mean particle spacing R s 
varies with v as (2ir/(v — l)) 1 ^ 2 assuming the square lattice structure. Now we construct 
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the phase diagram in the filling factor-temperature plane. The initial configurations are 
chosen randomly. We let the MD simulations evolve with time, make sure that the system 
reaches the equilibrium state, then freeze the skyrmion configuration and measure the static 
structure factor S(k) 

5( k ) = 1 ^ e *-CH*W-Ri(t)). (12 ) 

By measuring 5(k) and visually checking the snap shots, we identify the corresponding 
phases. First we have studied the zero temperature phase diagram. Insets of Fig. (2) demon- 
strate S'(k) at v — 1.02 and 1.07 exhibiting the six and four-fold peaks, as expected for the 
triangular and square lattice structures, respectively. At zero temperature, with the increase 
of filling factor away from v = 1, the system initially forms a triangular lattice and above a 
critical value of v\ = 1.025 a square lattice. At higher densities, the square Skyrme crystal 



will melt due to strong quantum fluctuations. [14,20] Subsequently, we investigate the melt 



ing lines for the various Skyrme crystals via finite temperature MD simulations. The melting 
phase boundary can be established by measuring the peak heights of static structure factor 
S'(Ko) as a function of temperature and K = (2tt, 2tt) for the square Skyrme crystal. The 
Lindemann criterion says that the solid melts when the thermal fluctuations of the equilib- 
rium positions are about 15% of the mean particle spacing R s , that is, V < u 2 > = clR s 
with cl — 0.15. Assuming Gaussian fluctuations, we obtain that S'(Ko) = e -K » <u2> / 2 = 1/3. 



21| Thus when the peak height of S'(Ko) reaches about 1/3, we consider the crystal having 
melted. In Fig.(U), we have plotted the melting phase boundary and noticed that at around 
v = 1.1, the melting temperature is about 0.035p s . The experimental value of p s is about 
1.5K yielding T m = 50mK, which is very close to the experimentally suggestive melting 
temperature. The spin structure factor (k) has also been measured 

S^(k) = — re^^W-^Wle'^W-^W). (13) 

Within the square Skyrme crystal phase, our simulations clearly demonstrate that the system 
possesses an antiferromagnetic spin order by exhibiting the peaks of S^(k) at k = (ir, n). |2(| 



We note that although it is possible that there could be a regime where the triangular and 
square lattices could compete to produce a disordered phase, we do not observe this with 
our model using the experimentally relevant parameters. In order to clarify this interesting 
possibility, further studies including careful finite size scaling analysis is required. [^(J In 
this work, we do not treat the terms in the action that guarantee the fermionic nature of 
skyrmions, nor do we include any quantum dynamics in the MD simulation. |22| Finally, 
we ignore any three-body interactions, treating the spin texture action as a set of pairwise 
interactions between skyrmions. 

In summary, we have determined an approximate two body interaction potential for 
skyrmions from an effective action via a variational calculation. Using this potential we 
have performed molecular dynamics simulations that have determined the phase diagram of 
the skryme system without having to posit a crystal structure. We find that the melting 
temperature within our model is close to the specific heat peaks seen in proposed skyrme 
crystal experiments, which gives a strong support to the very existence of skyrme crystal. 
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FIGURES 

FIG. 1. The two-body interaction potential as a function of R and Oq, as determined by vari- 
ational calculations of the total energy The curve with Oq = represents Vq (R) . Note the steep 
dependence of the energy as a function of twist. 

FIG. 2. Thermodynamic phase diagram of a manybody skyrmion system: Below v\ ~ 1.025, 
triangular lattice is stable. Above v\, square lattice is formed, which melts at around 0.035p s for 
v = 1.1. At some higher density, square skyrme crystal will melt due to quantum fluctuations. 
Insets show static structure factor S(k): At filling factor v = 1.02, six-fold peaks appear at the 
corresponding points to the triangular lattice. At filling factor v = 1.07, four- fold peaks appear at 
the corresponding points to the square lattice. 
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